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Abstract. The spin 1/2 Heisenberg model on a square lattice with antiferromagnetic nearest- and next-nearest neigh- 
bour interactions (the 71-/2 model) has long been studied as a paradigm of a two-dimensional frustrated quantum 
magnet. Only very recently, however, have the first experimental realisations of such systems been synthesized. The 
newest material, Pb2VO(P04)2 seems to have mixed ferro- and antiferromagnetic exchange couplings. In the light 
of this, we extend the semiclassical treatment of the J1-J2 model to include ferromagnetic interactions, and present an 
analysis of the finite temperature properties of the model based on the exact diagonalization of 8, 16 and 20 site clusters. 
We propose that diffuse neutron scattering can be used to resolve the ambiguity inherent in determining the ratio and 
sign of 71 and J2 from thermodynamic properties alone, and use a finite temperature Lanczos algorithm to make predic- 
tions for the relevant high temperature spin-spin correlation functions. The possibility of a spin-liquid phase occurring 
for ferromagnetic Ji is also briefly discussed. 

PACS. 71.27.-l-a Strongly correlated electron systems; heavy fermions - 71.10.-w Theory and models of many- 
electron systems - 75.40.Cx Static properties (order parameter, static susceptibility, heat capacities, critical exponents, 
etc.) 



1 Introduction 

The antiferromagnetic (AF) Heisenberg model on a square lat- 
tice with added nearest neighbour bonds, commonly refeiTed to 
as the Ji -J2 model, has long served as a paradigm for a two di- 
mensional frustrated AF. For a model with only one adjustable 
parameter - the ratio of next-nearest to nearest neighbour ex- 
change J2 jJi - it has an extremely rich phase diagram, offering 
the chance to study Neel order (NAF) with a reduced sublattice 
magnetization for J2l-I\ ^ 1, a collinear AF (CAF) phase se- 
lected by an order from disorder effect for 72 /-/i ^ 1, and a 
spin-gapped phase (or family of phases) for intermediate cou- 
pling. A correspondingly rich and occasionally controversial 
literature, has accompanied the development of these ideas 1 1 ] . 

More recently, the discovery of high superconductivity 
in doped layered cuprates whose undoped parent compounds 
are spin-half antiferromagnets, has lead to a renaissance of in- 
terest in 2D frustrated magnets. The 7i-72 model in particular 
has attracted renewed attention, both because of its simplicity 
and the possibility that its spin-gapped phase might provide a 
realization of Anderson's resonating valence bond (RVB) con- 
cept (see e.g. f2!|). 

Given all of this theoretical activity, it is surprising that the 
behaviour of the 7i -72 model in the presence of ferromagnetic 
(FM) exchange remains largely unexplored. 

Frustrated FM's have an interesting history in their own 
right. The solid phases of He III have very complex magnetic 



behaviour determined by competing FM and AF multiple spin 
exchange processes. Under appropriate conditions He III may 
exhibit FM, AF or spin liquid ground states. Another, solid 
state, example is provided by the doped colossal magnetoresis- 
tance (CMR) manganites, where the competition between ki- 
netic energy driven FM and superexchange driven AF is widely 
believed to result in a phase separation into regions with dif- 
ferent magnetic order. It is therefore worth asking whether the 
7i-72 model with FM 7i shows similarly exotic behaviour 

It is also perhaps a little surprising, considering the wide 
variety of magnetic materials now under study, that the first 
"7i-72 cornpound", Li2VOSi04, was discovered only very re- 
cently OB El- Li2VOSi04 is an insulating Vanadium oxide, 
with spin 5=1/2 ions arranged in square lattice planes, at 
the centres of VO4 pyramids. These are linked by Si04 tetrahe- 
dra, with Li ions occupying the space between the V-O planes. 
Because of this relatively complex structure, the magnetic ions 
are well separated, with weak superexchange between nearest 
and next-nearest neighbour V ions mediated by more than one 
intermediate O ion. A small number of related materials have 
now been synthesized, and the preliminary analysis of one of 
these, Pb2VO(P04)2 seems to offer evidence of mixed FM and 
AF exchange coupling |6]. However, a precise and unambigu- 
ous measurement of the exchange couplings 7i and 72 by, e. g., 
inelastic neutron scattering, has yet to be accomplished for any 
of these new compounds. 
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Motivated by Pb2VO(P04)2, in this paper we provide an 
overview of the ground state and finite temperature properties 
of the JrJ2 with mixed FM and AF coupHngs. We present a 
comprehensive semiclassical analysis of the three dominant or- 
dered phases of the model - a uniform FM phase, and q = 
{n,n) Neel (NAF) and q* = (7i;,0),(0,7i:) collinear (CAF) an- 
tiferromagnetic phases - together with an exact analytic diago- 
nalization of an eight site cluster, and finite temperature Lanc- 
zos results for the heat capacity and magnetic susceptibility for 
16 and 20 site clusters. 

We argue that, in addition to the known spin liquid region 
for Ji > 0, J2 ^ Ji /2, where the NAF and CAF phases com- 
pete, a new spin liquid region may exist for Ji < 0, ^2 ~ /2, 
where the FM and CAF phases compete. We also propose that, 
because of their low magnetic energy scales, diffuse neutron 
scattering at finite temperatures can provide a very useful source 
of information about the nature of the competing magnetic in- 
teractions in these materials. With this in mind, we present the 
first quantitative numerical estimates of the magnetic structure 
factor S(^, T) for the J1-J2 model. 

2 Zero temperature properties 

2.1 Classical phase diagram and general arguments 

We consider the spin 1/2 Heisenberg model on a square lattice 

J^^JiY.SrSj+J2Y.Si-Sk (1) 

{'■;■> 1 («>2 

where the sum on runs over nearest neighbour and the 
sum {ik)2 over diagonal next-nearest neighbour bonds. We al- 
low the exchange constants Ji and J2 to be negative (FM) as 
well as positive (AF) 

Since the properties of the Ji -J2 model depend on the rela- 
tive and not the absolute size of the exchange couplings Ji and 
J2, it is convenient to characterize it by an overall energy scale 

Jc - \/jf+Ji (2) 

and a frustration angle such that 

Ji ^ JcCos{(j)) J2~JcSin{^) 

=tan-'(72M) (3) 

As discussed in Section|3] 7c can in principle be determined di- 
rectly from the asymptotic behaviour of heat capacity and sus- 
ceptibility at high temperatures. However the different physical 
properties of the model depend on the angle (j), and this is much 
harder to determine from experiment. 

Let us first consider the simplest possible classical analysis 
of the model. We assume that the system orders at zero tem- 
perature in such a way that all the spins are oriented in a com- 
mon plane |7]. In this case the ground state energy of the J1-J2 
model is minimised by an order parameter with wave vector 
q = q* such that the energy per spin 

E{q*) = ]^zS^[Jiy,{q*)+J2y2{q*)] (4) 



J2 




J2/J1 = 



Fig. 1. Classical phase diagram. Numbers are ratios of exchange cou- 
plings J2/J\ for phase boundaries as determined from the classical 
ground state energy. The boundary between FM and NAF phase is 
the line J\ = J2 < Q {Ji/Ji = in the figure). Values of J2/J1 
in parentheses show where zero point fluctuations destroy the relevant 
order parameter at a semiclassical level, as discussed in Section 
The shaded areas for Ji > correspond to the known spin-liquid 
regime, and for 7] < to another spin liquid region. The frustration 
angle is given by (j) = tan^' (72/71 ). 
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Fig. 2. Classical energies E{q*) as a function of the frustration angle 
(j) in units of 7c. Solid line: FM, q* = 0; dotted line: NAF, q* = [tz, 7t); 
dashed line: CAF, q* = (;r,0); long-dashed line: four sublattice state 
for^* = (n:/2,0). 



takes on its minimal value. Here 

Yi{q) (cos(^;t)+cos(^_v))/2 

72(^) cos(^.^)cos(g'y) (5) 

z = 4 is the lattice coordination number for each type of bond 
and S = 1/2 is the size of the spin. 
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q* 


Energy 


Range (7i,72) 


Range ^ 


NAF 


{n,7i) 


-J1/2 + J2/2 


yi>0, J2<Ji/2 


-nl2 <(j) < tan"' (5) 


CAF 


(0,k) or (;r,0) 


-J2/2 


\J2\>\Jl\/2 


tan"' (5) <(j> < n-lnn^^{j) 


FM 


(0,0) 


J1/2 + J2/2 


yi<o, j2<-J\/2 


;r-tan"'(3) < ^ < -7t/2 



Table 1. Parameters for classical ground states diagram. 



This analysis selects three phases, a Neel AF (NAF), a collinear 
AF (CAF) and a uniform FM to give the phase diagram shown 
in Figure n with the parameters given in Tabled Note that the 
coplanar spiral states with q* = (2;r«/m,0), where {n,m} are 
integers, have energies which interpolate between the CAF and 
FM, and all of these states become degenerate exactly at the 
transition from FM to CAF. The classical energies of FM, NAF 
and CAF order parameters, together with a four sublattice state 
with q* — {7t/2,Q) are shown in Figure|2] 

Each of the three lines selected as classical phase bound- 
aries have interesting properties which are related to symme- 
tries of the Hamiltonian. The mirror symmetry of the classical 
phase diagram about the line = is particularly easy to un- 
derstand. Since the square lattice is bipartite, the spins can be 
divided into A and B sublattices and, the Hamiltonian remains 
invariant under the transformation 



(6) 



This converts the classical NAF into a FM, and the classical 
q* = (0,7r) CAF into the classical q* = (7r,0) CAF. Exactly 
on the line 7i = the A and B sublattices are entirely discon- 
nected, so it is possible to rotate the classical NAF into the FM, 
and the classical q* — {7t,0) CAF into the classical q* = (0, 7t) 
CAF, without any cost of energy. 

Of course it is reasonable to ask whether such a naive clas- 
sical picture has any relevance at all for the physics of a two 
dimensional frustrated spin 1/2 magnet. By way of an answer, 
let us consider in turn the simplest possible quantum analysis 
of the model. 

Using the geometric trick of double counting all J\ bonds 
and then setting Ji Ji /2, the cross-linked square lattice can 
be expressed as four interpenetrating sublattices of edge shar- 
ing tetrahedra (cross-linked squares). And since each of these 
cross-linked squares is a complete graph, the Hamiltonian 
can be rewritten in terms of the total spin on a square 



Jf ^-2/28(8 +l)N- 



□ 



□ 



(7) 



where the sum runs over all squares q of the lattice, and 

[£2j 



4 + 



h 



h 
2 



(8) 



The spins within each square are counted clockwise (or coun- 
terclockwise) such that 



X2b = 82 +8^ 



Si 



(9) 
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Fig. 3. On the transition lines a) J2/J] = 0.5 and b) J2/Ji = —0.5 
the Ising domain wall energy between CAF and NAF or FM states 
vanishes and the system can break up into stripe like domains with all 
possible values of total magnetization. 



With the Hamiltonian written in this way, the special role of the 
classical transition lines J2±Ji/2 — becomes self-evident. 
Moreover, on these lines, the lowest energy is achieved by any 
state for which the sum (difference) of the A and B sublattice 
magnetisations vanishes in each square. 

Classically, this condition can be satisfied by an exponen- 
tially large number of states. As a result, at these transitions, 
it becomes possible to break the system up into stripes. In the 
Ising limit these stripes are the alternating domains illustrated 
in Figure |3] The special (local) degeneracy of the Hamilto- 
nian (0 also reveals itself in lines of zeros in the spinwave 
spectrum, discussed below. 

We now use equation (0 to construct a minimal quantum 
theory for the model. The state of the whole system is deter- 
mined if we specify all of the spins on any of the four sublat- 
tices of tetrahedra. So a lower quantum bound on the ground 
state energy can be obtained by considering a single isolated 
tetrahedron. 

We find what are essentially the same three ground states 
as a function of 0, with energies and quantum numbers given 
in Table 12] The full set of energy eigenvalues is shown in Fig- 
ure|3 Apart from in the FM phase, where they must agree, the 
upper bound on the ground state energy obtained from the sim- 
ple classical analysis, and the lower bound obtained from this 
minimal quantum estimate are generally quite different. How- 
ever, it is far more significant that the same three phases are 
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Energy 


Range (Ji ,^2) 


Range (p 


"NAF" 


(0,1,1) 


-A +J2/2 


71 >0, J2<Ji/2 


~7t/2 <(j)< tan-^ij) 


"CAF" 


(0,0,0) 


-372/2 


J2>Jl/2 


tan"' (2) <(j> < n-tan'^d) 


FM 


(2,1,1) 


J1/2 + J2/2 


Ji<0, J2<-Ji/4 


;r-tan"'(|) < (j) < -n/2 



Table 2. Parameters for the ground states of a single cross linked square. 
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Fig. 4. Energies of the six eigenstates of a single tetrahedron, as a 
function of the frustration angle (j), in units of Jc- 



found', with the same phase boundaries - excepting that be- 
tween the CAF and FM, where quantum fluctuations extend 
the CAF regime at the expense of the FM regime. 

This simple correspondence between the most naive clas- 
sical and quantum theories lends us some confidence in both. 
And although neither can give a complete description of the 
model, both can be improved systematically, as discussed in 
later sections of the paper The key question which remains is 
how the phase transitions between the three dominant phases 
take place. The existence of an extensive set of spiral states 
degenerate with the FM and CAF order parameters at (j) = 
TT — tan"' (1/2), and the complicated level crossings for the 
tetrahedron near (p = tan"' (1 /2) and (p = n — tan"' (1 /2) al- 
ready hint that these can be non-trivial. 



2.2 Semiclassical spin wave analysis 

The first step to improving on the naive classical phase dia- 
gram is to consider the influence of semiclassical spin wave 
excitations. While for the FM, spin waves are eigenstates, the 
frustrated NAF and CAF phases both show marked fluctuations 
at zero temperature. All three phases must be unstable at finite 

' Here "phase" should be understood to mean a ground state wave 
function with the same total spin and spatial symmetries as the corre- 
sponding classical order parameter — SU(2) invariance clearly cannot 
be spontaneously broken in a small cluster. 
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Fig. 5. Evolution of the spin stiffness D in the FM phase (solid line), 
and of the spinwave velocities v in the NAF phase (dotted line) and v^, 
Vy in the CAF phases (dashed lines), as a function of the frustration 
angle ^ in units such that the lattice constant a = 1 and J^ = \. 



temperatures, as long range order would violate the Mermin- 
Wagner theorem. None the less, linear spin wave theory cap- 
tures the essential physics of the ordered phases, and provides 
some interesting hints about how the classical phase diagram 
must be modified in the quantum case. 



2.2.1 FM phase 

Expanding about the FM phase we find a spin wave dispersion 



(S){ci) 
where 



-45 [/i + 72] + 45 [7i (c, + Cy) /2 + hcxCy\ ( 1 2) 



Cv = COs(^;t) Cy = COs(g'v) 



(13) 



in units such that the lattice constant a ~ I. The spin wave dis- 
persion for a range of values of throughout the FM phase is 
shown in Figure|6l Note that in this and all subsequent plots of 
spin wave dispersion, the and qy values run from — tt to n, 
i. e. over the full Brillouin zone (BZ) for the square lattice, and 
not the reduced magnetic BZ's appropriate to the NAF or CAF 
phases. 

While the fully polarized FM ground state remains an exact 
eigenstate of the frustrated model, its dispersion is profoundly 
modified by competing interactions. At the boundary with the 
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Fig. 6. Evolution of spinwave dispersion in FM phase. From left to right 
exchange, border with CAE 



border with NAF, deep within FM phase, pure nearest neighbour 
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Fig. 7. Evolution of the heat capacity Cy as a function of the frus- 
tration angle (j). In FM regions the quantity plotted is the prefactor y 
of Cv = yT, and in AF regions the prefactor /3 of Cv = jST^, where 
temperature is measured in units of J^. 



NAF phase for (pc = —7^/2, the dispersion is that of a pure J2 
FM, which has the same magnetic BZ as the NAF phase, and 
therefore zeroes at ^ = (tt, n) and symmetry points, in addition 
the usual <7 = (0, 0) Goldstone mode. Deep within the FM phase 
for (j) — — ;r + tan^' (1 /2), the dispersion behaves as 



CO{q) ^ Dq^ 
where the stiffness constant D is given by 

D = -{Ji+2J2)S 



(14) 



Cv = 



C(2) 
2k \D 



(16) 



where ^(2) = The coefficient of T as function of is 

plotted in Figure It diverges at the transition between the 
FM and the CAF, but approaches a constant at the transition 
between FM and NAF. 



2.2.2 NAF phase 

The spinwave spectrum for the NAF is given by 
CO{q) 



(17) 



where the coupling between spins on a given sublattice is 

A,-45[yi-y2(l-c,Cv)] (18) 
and the coupling between the two sublattices is 

B„ = 2/1 5(c,- + c,) (19) 

Where J2 is FM, it acts to stabilize the NAF order, where J2 is 
AF, it acts to destroy it. Once again this competition is visible 
in the spin wave dispersion, as show in Figure|8] 

At the boundary with the FM phase for ~ —7^/2, the 
dispersion is that of a pure J2 NAF, and exactly matches that 
of the FM on this phase boundary. Deep within the NAF phase 
for (pc = — tan^'(l/2), the low energy spin wave dispersion 
behaves as 

COiq)^v,\q-q*\ (20) 



(15) where the isotropic spin wave velocity is given by 



in the zone centre, but vanishes on the zone boundary. The vari- 
ation of D as a function of (j) is plotted in Figure|5lFor (j) — — n 
the dispersion is that of the familiar pure Ji FM. And, finally, on 
the boundary between FM and CAF for 0^ = TT — tan^^ (1 /2), 
the dispersion vanishes on the lines qx — and qy — 0. These 
lines of zeros are a direct manifestation of the special local 
symmetry discussed in Section lzTI 

The heat capacity of a FM in 2D is linear at low tempera- 
tures, reflecting a constant density of states at zero energy, and 



v, = 2S^2Ji{Ji~2J2) (21) 

and q* ~ {7Z,7z), as expected. However it exhibits a marked dis- 
persion about the magnetic zone boundary, as compared to the 
pure Ji NAF for ^ = 0. Finally, on the boundary with the CAF 
phase for (j)c = tan^^ (1 /2), the dispersion vanishes on the lines 
qx = 0, qx = ±7t and q^. ~ 0, q^ ^ ±n. Values of the spinwave 
velocity V5 are shown in Figure |5j in units in such that the lat- 
tice spacing a = 1 and the overall energy scale 7c = 1 • 
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Fig. 8. Evolution of spinwave dispersion in NAF phase. From left to right - 
exchange, border with CAE 



border with FM, deep within NAF phase, pure nearest neighbour 




Fig. 9. Different estimates of the ground state energy per spin as a 
function of the frustration angle (j) in units of Jc. Uppermost (solid) 
line - classical energy of FM, NAF and CAF order parameters given 
by Equation|4] Lowermost (dotted) line - lower quantum bound given 
by the ground state of the Hamiltonian 0. Inner (dashed) line - semi- 
classical spinwave theory. Diamonds - variational lower bound based 
on diagonalization of a nine-site cluster. Solid and open circles - 
ground state energy obtained from exact diagonalization of a 16- and 
20-site cluster, respectively. Note that in the collinear phase the vari- 
ational bound from the nine-site cluster and the exact-diagonalization 
result for the 16-site cluster nearly coincide. 



Fig. 10. Evolution of the inverse transverse susceptibility a 
function of the frustration angle (p, in units such that Jc = I and 



1. 



estimates It also lies very close to our numerical estimates 
and an improved variational quantum bound, discussed below. 
These are also plotted in Figure|9l 

It is also interesting to consider the dependence of the 
experimentally accessible transverse susceptibility X± which, 
within spin wave theory for a two sublattice AF is given by 



The low temperature heat capacity of the NAF is controlled 
by the spin wave velocity and is given by 



Cv 



3C(3) 

27Z 



(22) 



where ^(3) = 1.202.... The relevant coefficient of T^, as a 
function of is plotted in Figure0 

The quantum zero point corrections to the ground state en- 
ergy also vary strongly as a function of (j) within the NAF 
state, vanishing altogether at the boundary with the FM. These 
are shown in Figure |9] together with the upper classical and 
lower quantum bounds discussed in Section IZD In fact, near 
the transition from NAF to CAF, the semiclassical estimate of 
the ground state energy surprisingly well with series expansion 



(23) 



in units such that the overall prefactor {g^sY ~ 1- The varia- 
tion of the inverse susceptibility as a function of ^, normalized 
to the value for = 0, is shown in FigurefTol The susceptibiUty 
is a continuous function of 0, diverging (as Xx ~ (0 — 0*)^') 
at the FM phase boundary, and matching that of the CAF at the 
other end of the NAF phase. However its derivative is discon- 
tinuous across each transition, reflecting the changing symme- 
try of the order parameter. 

We can gain still more information about the evolution of 
the NAF state for different couplings by calculating the sublat- 
tice magnetisation m^. Quantum fluctuations reduce from 
its classical value Ma — 1/2, and in terms of the coefficients 
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Fig. 11. Evolution of sublattice magnetization iiis as a function of the 
frustration angle <j). 




defined above, it is given by: 



(24) 



The variation of as a function of (j) is shown in Figure fTTI 
The sublattice magnetization of the NAF vanishes for (j)/n = 
0.12, corresponding to a ratio of Jo/Ji = 0.38. This result im- 
plies that, at a semiclassical level, the NAF order parameter is 
destroyed by fluctuations long before the competing CAF be- 
comes energetically favourable. Historically, this was the first 
signature of the existence of an intermediate, spin-gapped phase 
between the NAF and CAF in the purely AF Ji -J2 model 1 9 ] . 
We return to this point in the context of the model with mixed 
FM and AF couplings below. 



2.2.3 CAF Phase 

The CAF phase can be thought of as two interpenetrating NAF 
lattices of diagonal J2 bonds. At a purely classical level, these 
sublattices are decoupled and can be rotated freely about one 
another. However quantum fluctuations stabilise the configura- 
tion in which the Neel vectors of both sublattices, and therefore 
the associated spins, are collinear. This tendency of fluctuations 
to favour collinear spin configurations is well known 1 10, 11, 
[T2il . and for small J1/J2 is independant of the sign of Ji. We 
have checked explicitly that this order from disorder effect sur- 
vives for larger, FM values of Ji by peforming spin wave calcu- 
lations for more general four-sublattice states. Details of these 
will be reproduced elsewhere. We note that any further two- 
sublattice canting of the CAF state is energeticaly unfavourable 
at a classical level, even for quite large FM Ji . 

Once the two sublattice CAF has been selected, the analy- 
sis of the spinwave spectrum is straightforward. The spinwave 
spectrum for the CAF is once again of the form equation 



For the CAF order parameter with ^* (tt, 0) we find 

A,^2S[2J2+JiCy] (25) 
Bg =2Sc,[Ji+2J2Cy] (26) 

in accordance with lH (die result for q* ^ (0, n) can be ob- 
tained simply by exchanging x and y above). The evolution of 
the spin wave dispersion within the CAF phase is shown in 
Figure|S] plotted within the full BZ for the square lattice. 

At the border with the NAF phase, for (j)^. ~ tan"^ (1 /2) the 
CAF has lines of zero modes for = 0, q^ ~ ±7r and qy — ±7t, 
but not for qy = 0, and has maxima at ^ = (±7r/2,0). Within 
the CAF phase for AF couplings, for (p ~ 7t/4, the dispersion 
has peaked maxima at ^ = (±;r/2,0), dispersionless ridges for 
qy — ±n/2. In all cases the dispersion had zeros at the wave 
vectors appropriate to the order parameter, i.e. q = (0,0) and 
q* = (±;r,0). Near to these we find a linear but anisotropic 
spin wave dispersion 



where the spin wave velocities are given by 



2J2+J1 



(27) 

(28) 
(29) 



= ^i2J2+Ji)i2J2-Ji) 

The low temperature heat capacity of the CAF is controlled by 
the average spin wave velo 
equivalent of equation il6\ 

3C(3) 



the average spin wave velocity v'j=(vivv)i and is given by the 



Cv 



2% 



TV 
vs) 



(30) 



where ^(3) = 1.202.... The relevant coefficient of T^, as a 
function of <j) is plotted in Figure 

The symmetry between the x and y axes is broken by the 
CAF order parameter, which is reflected in the different values 
of Vv and Vy. However, at a semiclassical level, this symme- 
try breaking is not reflected in different values the transverse 
susceptibility x±x and x±y - plotted in Figure^] This can be 
understood as follows - the transverse susceptibility associated 
with an AF can be expressed in terms of the spin stiffness Ps 
and the spinwave velocity using the hydrodynamic relation 
X± = Pi/i^'.?- However at this level of approximation the varia- 
tion of D with angle in the plane is precisely that required to 
cancel the variation of v^. Once again, the special role of the 
lines Ji = ±2^2 is evident - at the borders of the CAF phase 
the solutions for Vy become imaginary, while the transverse sus- 
ceptibility diverges. 

For the pure next-nearest neighbour model at (j) = n/2, 
the magnetic BZ is further reduced, with additional minima 
at ^ = (±7r,±7i;) and qy — ±n. The dispersion exhibits ridge 
like maxima for q — {±n /2,zLn /2). Within the CAF phase for 
partially FM coupling at (j) ~ 37z/4, the dispersion is that for 
(j) — 7z/4 described above, but with q translated by ( ;r/2 , 0) . Fi- 
nally, for 0c = ?r — tan^'(l /2), at the boundary with the FM 
state, it has maxima for q = (±;r/2,±;r) and line zeros for 
qx = ±7t and qy ~ 0, qy = ±n, but not for q^ = 0. 

While it was possible to exactly match up the spin wave 
dispersion on the boundary between the NAF and FM, where 
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the A and B sublattices decouple, this cannot be done for the 
boundaries between the CAF and NAF or CAF and FM. The 
different way in which the CAF order parameter breaks the lat- 
tice symmetry is immediately apparent in the different number 
of lines of zero modes present in the CAF for (j) — tan^' (1 /2) 
and ^ = K — tan^' (1 /2), as compared with the NAF or FM. 

This rules out a smooth transition from one state to the 
other, and at a semiclassical level, the only way in which the 
system can avoid embarrassment is to dissolve the classical 
CAF order parameter before the phase boundary is reached. 
And indeed the relevant sublattice magnetisation does vanish 
for (p/n = 0.15 (J2/J1 = 0.52), as the CAF approaches the 
NAF, and for (p/7Z = 0.85 (J2/J1 = -0.51), as the CAF ap- 
proaches the FM phase. Consistent with this, the heat capacity 
(Figure shows clear anomalies as a function of <j) at either 
end of the CAF phase. 



2.3 Beyond the semiclassical picture 

The simple quantum and semiclassical arguments presented 
above give rise to an equally simple and self consistent pic- 
ture of the phase diagram of the J1-J2 model as a function of 
(j). The model has three dominant phases, FM, NAF and CAF. 
The phase transition between the FM and NAF is straightfor- 
ward. The phase transitions between the CAF and NAF, and 
the CAF and FM are not, and probably, take place through an 
intermediate phase. 

The different estimates of the ground state energy of the 
system in Figure |9] give us further reason to believe that this 
semiclassical picture is not far from the truth. Clearly the spin 
wave estimate of the ground state energy lies a long way above 
the lower bound from the tetrahedral cluster. But, as discussed 
in Appendix |C] it is possible to construct a much better quan- 
tum bound variationally, from the numerical diagonalization of 
a 9 site clusters with modified boundary conditions. Increas- 
ing the cluster size to 13 sites does not lead to any significant 
change in the ground state energy. Since the semiclassical es- 
timate is within a few percent of the quantum bound on the 
ground state energy for almost all values of (j), there is little 
room for drastic changes in the phase diagram. 

But of course this is not the end of the story. Even if our 
semiclassical estimate of zero point energy is reliable, fluc- 
tuations will strongly renormalise the semiclassical spin wave 
spectra and correlation functions. In the pure Ji Neel AF L13il 
these coiTections lead to an enhancement of the spin wave ve- 
locity 



Vs ZcVs Zc = l 



0.1580 0.0216 



2S 



(25) 



2 



1.1794 



a suppression of the perpendicular susceptibility 

, 0.551 0.065 

X,^Z,x. Z, = l-^ + ^ + ...«0.514 

and a slight enhancement of the sublattice magnetisation. 



ms 



' ZjjiS 



7 — 1 — ■ 



0.383 0.007 



25 



(25)2 



! 0.613 



The values of the coefficients Z as a function of Ji and J2 for 
^2 7^ remain to be calculated, and until they are known, quan- 
titative comparison of the numbers found above with experi- 
ment must be approached with some caution. However experi- 
ence with the pure Ji model suggests that, as long as fluctua- 
tions do not destroy magnetic order completely, coiTections are 
reasonably uniform and the semiclassical description of spin 
correlations remains qualitatively, if not quantitatively valid. 

This "renormalised classical" physics should be expected to 
break down near the highly frustrated transitions into, and out 
of the CAF phase. The large body of existing work on the J1-J2 
model with AF couplings suggests that the space between the 
NAF and CAF phases is filled by a spin-gapped 'spin-liquid' 
phase. Spin liquids are known to occur adjacent to a FM phase 
in the Heisenberg model on a triangular lattice with competing 
FM and AF cyclic exchanges 1 14], and on a honey-comb lat- 
tice with competing Ji , J2 and (next-next nearest neighbour) 
interactions 1 15]. 

So what happens between the FM and CAF phases in our 
model ? 

At least at first glance, the situation seems to be very simi- 
lar A high local degeneracy in the classical spectrum - a family 
of degenerate order parameters for (p = n — tan^ ' (1 /2) - leads 
to line zeros in the spin wave velocity and vanishing sublattice 
magnetisation exiting the CAF state towards the FM. Exactly 
the same things happen at the much studied boundary with the 
NAF. This similarity is by no means proof of the existence of 
a spin liquid state in the J1-J2 model with FM 7i, but it is a 
sufficient reason to start looking for one. 

Further circumstantial evidence in favour of this hypothe- 
sis can be obtained from the exact analytic diagonalization of 
an 8-site cluster described in Appendix 151 The resulting en- 
ergy spectrum, classified by spin, is shown in Figure [O] The 
straight-forward phase transition between NAF and FM states 
for (j) — — n/2 shows up as a multiple crossings of ground state 
and excitation energy levels, all of which take place at the same 
critical value of = — ;r/2. 

Where the singlets associated with NAF and CAF order pa- 
rameters cross, the reordering of excited states does not take 
place at a single critical value of <p, but is spread out over a finite 
range of 0. Simply counting where the lowest lying triplet ex- 
citation crosses the lowest lying singlet excitation either side of 
the ground state crossing gives a remarkably good (if arbitrary) 
estimate of the extent of the spin liquid region - from Ji /J2 = 
0.38 to J1/J2 = 0.60, values which are comparable with those 
found in the existing literature ]2, 5..9> .ld..l7..1&..19t.2Qil . 

Examining the level crossings associated with the transi- 
tion from CAF to FM we see the same extended structure. In 
this case applying the same naive criterion based on the cross- 
ing of first excitations would predict a spin liquid region from 
J1/J2 = —0.38 to J1/J2 = —0.60. However, in order to obtain 
a more serious numerical estimate of the domain of stability 
of the CAF and FM order parameters it would be necessary to 
look at the finite size scaling of not only the first excitation en- 
ergies, but of the entire "Anderson Tower" of states i2lll23l23ll 
which go to make up the CAF order parameter in the thermody- 
namic limit, for a sequence of clusters including those of rela- 
tively large sizes (e. g. 32 and 36 sites). This analysis is beyond 
the scope of the present paper. 
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Fig. 13. Energy levels (per spin) of the 8 site cluster, classified according to total spin 12 as a function of the frustration angle (j), in units of J^. 
solid lines ~ Q = (singlet); dotted lines - £2 = I (triplet); short-dashed lines - Q = 2; long-dashed lines - Q = 3; dash-dotted line - Q = 4 
(maximal spin). 



It also seems premature to speculate about the nature of 
any possible new phase appearing between the FM and CAR 
Thinking classically, one might imagine that a spiral or canted 
state arises which interpolates between the CAF and FM. Adding 
an AF Jj, interaction to the model would tend to favour such 
states (cf. (T^). However our preliminary analysis suggests 
that, for the relevant range of Ji and J2 (with = 0), these 
states are still more unstable against fluctuations than the CAF 
with which they compete. 

In short, while the outcome remains uncertain, there is clearly 
reason to suspect that something interesting happens at the tran- 
sition from CAF to FM. In marked contrast to the J1-J2 model 
with AF interactions, the exis ting literature on this problem ap- 
pears to be in its infancy 12411251 . 



3 Finite temperature properties 

3.1 General considerations 

The energy and temperature scales of the Ji -J2 model are con- 
trolled by the single parameter J^. In principle, this can be de- 
termined directly from the asymptotic behaviour of the mag- 
netic contribution to the heat capacity at high temperatures 113. 



CviT- 



3Jl_ 

8 r2 



(31) 



written here in units 'natural' such that = 1 . It also controls 
the deviation of the high temperature magnetic susceptibility 
from a simple Curie law 



;^-'(r->oo) 



- ©cw + 5 



c 



(32) 



where we again work in 'natural' units such that {gUs)^ = 1 
S{S+l) 



c 



1 

4 



and 



0cw=4(7i+72)C = /i+/2. 



(33) 



(34) 



However, while knowledge of both Jc and ©cw fixes Ji +J2 
and I 7i — ^2 I . the sign of Ji — J2 remains undetermined since 
there are two possible values (j)± of the angle (j). And because (j) 
determines the physics of the J1-J2 model, this uncertainty can 
lead to alternative parameterizations of the model which lie in 
completely different phases. 
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Pb2VO(P04)2 




Li2VOSi04 




Li2VOGe04 


















lil 


LL8] 


0cw[K] 


4 




7.4 


8.2 


9.65 


7.2 




5.2 


9.8 


0cw/7i 


0.49 




1.39 


1.69 








1.49 






0.67 


0.64 


0.41 


0.27 


0.47 


0.43 


0.36 


0.38 


0.43 0.33 






-0.11 


0.03 


0.13 






0.06 


0.08 


0.07 


(Ji/Jx)- 


-1.64 


-2.21 


3.5 


1.1 


11.7 


4.76 


2.13 


2.5 


4.76 1.69 


{J2/Jl) + 




-0.37 


0.1 


0.44 






0.18 


0.25 


0.24 



Table 3. Compilation of the experimental results and theoretical estimates on the Curie-Weiss temperature ©cw ~ {Ji +J2)/kB, the ratio 
©Cw/T/ of it to the maximum position of the uniform magnetic susceptibility x{T) the corresponding frustration parameters. The data are 
taken from 0HJg]. The displayed theoretical values for the frustration parameters obtained by fits to high-temperature series expansions are 
taken from llMIl^ . The reference numbers are used to label the corresponding columns. The unlabelled columns contain our own estimates 
derived from the dependence of &cy^ /T^ on 0, see Figure fTsl The ± subscripts of ^ and J2/J\ distinguish the two different possible points in 
the (^1,^2) phase diagram. 



The coefficients of the high temperature series expansions 
of the heat capacity and susceptibility J\ -J2 model are known to 
high order [18, 26], and their reliability at low temperatures can 
be greatly improved by carefully constraining the analytic con- 
tinuation of the high temperature series 1 19, 27]. None the less 
it proves very difficult to determine J\ and J2 unambiguously 
from experimental measurements of heat capacity and mag- 
netic susceptibility. In the case of Li2VOSi04, where both J\ 
and J2 are believed to be antiferromagnetic, estimates of the ra- 
tio of J2/J\ vary by more than a factor ten ifllHiElil-llll, see 
Table|3l Preliminary analysis of Pb2 VO(P04)2, where either J\ 
or J2 is believed to be ferromagnetic, does not unambiguously 
determine which is the ferro- and which the antiferromagnetic 
coupling 1 6]. 

In the light of this uncertainty and controversy, we have 
used our analytic solution of the 8-site cluster and an imple- 
mentation of the finite temperature Lanczos method (FTLM) 
to determine the heat capacity, uniform magnetic susceptibil- 
ity and static spin structure factor for the model on 16- and 
20-site clusters with periodic boundary conditions. In particu- 
lar, below, we present predictions for diffuse neutron scattering 
cross sections which are difficult to access by series expansion, 
and which can in principle remove all ambiguity about the pa- 
rameterization of the model. 

Since these calculations are based on small clusters of spins, 
they are of limited use in addressing questions such as the finite 
temperature Ising transition observed in Monte Carlo simula- 
tions of the classical square-lattice 71-/2 model |28]. However 
they should provide a reliable guide to physics at finite tem- 
peratures and short length scales, particularly in the frustrated 
phases of the model, where long range order is greatly sup- 
pressed. 

The FTLM is based on the Lanczos algorithm which is used 
to iteratively exactly diagonalize the Hamiltonian matrix for 
the cluster considered: The Lanczos algorithm starts with a 
randomly chosen vector in the Hilbert space. Successive ap- 
plications of the Hamiltonian eventually "rotates" the starting 
vector to the ground state of the system (if the starting vector 
is not orthogonal to it), thereby generating a tridiagonal matrix 
having eigenvalues and eigenvectors corresponding to the low- 
est eigenvalues and eigenvectors of the full Hamiltonian (this 
always has to be checked). The FTLM utilizes the eigenvalues 
and eigenvectors of 1000 successive diagonahzations with dif- 



000000000000 




Fig. 14. Tiles of size eight, 16, and 20 used in the finite-temperature 
calculations. In the lower right corner, the labelling of the two ex- 
change constants is illustrated. 



ferent random starting vectors to generate the partition function 
and the expectation values as defined below. For a detailed de- 
scription of it, see i29ll . 



The three clusters, together with the realisation of the spin 
exchange interactions, are shown in Figure^^ Due to the sym- 
metry of the model, these are the only possible squares of size 
with 4 < < 32 which are compatible with both coUinear 
and Neel order With three cluster sizes at hand, knowing that 
the eight-site cluster is almost a complete graph and the 16- 
site cluster a four-dimensional hypercube, we did not attempt 
to perform a finite-size scaling analysis. For this to be mean- 
ingful, results for a system size at least = 32, and as well as 
for different (i. e. open) boundary conditions would be needed. 
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Fig. 15. Ratio of the Curie- Weiss temperature ©cw to the maximum 
position of the uniform magnetic susceptibility x{T) as a function 
of the frustration angle (j). The solid (dashed) line denote the results 
for the 20-site (16-site) cluster, the dotted line comprises the result 
for the eight-site cluster The straight horizontal lines correspond to 
the experimental values ©cw/^i = 0.49, 1.39, 1.49 for Pb2VO(P04)2 
(dash-dotted) 01, Li2VOSi04 (dashed) 01, and Li2VOGe04 (long- 
dashed) PI, respectively. 



3.2 Heat capacity and magnetic susceptibility 

In units with dimensions restored, the heat capacity and the 
magnetic susceptibility are defined by 



XiT) 
Cv{T) 



Mb r2 



{H^)-{Hy 



(35) 



(36) 



where (. . .) denotes the thermal average, r^^'^-the z compo- 
nent of the total momentum of the system, and the number 
of sites of the system considered. Np^ is the Avogadro constant, 
/Xo the magnetic permeability, g the gyromagnetic ratio, /Xb the 
Bohr magneton, and the Boltzmann constant. For the present 
nonmagnetic (zero field) case we have (5'°') ~ 0. 

A commonly used experimentally accessible parameter to 
determine the frustration angle (and hence the value Ji/Jx) 
is the ratio of the position T^^ of the maximum of the magnetic 
susceptibility x{T) to the Curie-Weiss temperature ©cw- To 
avoid singularities, we plot the inverse of this quantity, which 
is shown in FigurefTsl 

Apart from the strongly frustrated spin-liquid region around 
^ « n/6, the differences in ©cw/T^ for the 16- and 20-site 
cluster are small. This suggests that the behaviour of x{T), at 
least for k^T > 7c, is dominated by correlations which are fully 
taken into account already by the small clusters, and therefore 
finite-size effects do not play an important role. 

A common feature of ©cw / 7^ for all cluster sizes is the ex- 
istence of two possible angles 0_ and 0+ for a given value for 
©cw/^x' reflecting the fact that from the knowledge of ©cw 
and 7c alone, cannot be determined unambiguously. We have 



indicated the experimental values taken from 0] and jal for 
©cw/7i for Pb2VO(P04)2, Li2VOSi04, and Li2VOGe04 by 
the thin horizontal lines in Figure^] For all three compounds, 
0- corresponds to a realization of the phase with strong collinear 
antiferromagnetic correlations, and 0+ to the phase where Neel- 
type correlations dominate. Interestingly, 
some values of 0+ lie very close to the spin-liquid regime 
0.115 <0/;/r< 0.183. 

Table |3l holds a summary of the values for ©cw / and 
0± = tan"' ((-/2M)±) found in the literature K,"", 18,131, 
together with our findings, ©cw = (-^i + Ji) /^b can be deter- 
mined by a fit of a Curie-Weiss law to the high-temperature 
tail of the susceptibility ;i^(r), as was done for the two Li com- 
pounds |4, 5]. This procedure can be improved by including 
higher-order terms in an expansion of x{T), see equation M>2\ . 
For the Pb compound |6], the high-temperature series expan- 
sion for found in |18] was applied. For the compound 
Li2VOSi04 fl|5(], this leads to errors of the order of 10 %. 

In contrast to ©cw, the determination of the frustration an- 
gle is much more involved. Melzi et al. 0<|1|] use exact diag- 
onalization data for the heat capacity of the eight- and 16-site 
clusters [30, 31], while again fits to high-temperature series ex- 
pansions are used in llfil fl3 fioll . In particular for Li2VOSi04, 
the results for 72/7i differ by more than an order of magnitude. 
Still they indicate qualitatively the same ordered phase in the 
ground state of the compound. 

However, 7i and 72 are not uniquely determined by this 
analysis of the temperature dependence of the magnetic sus- 
ceptibility or the comparison with the behaviour of the specific 
heat alone. We will return to this issue in the next section dis- 
cussing the static spin structure factor. 

We have computed the heat capacity Cy [T) in the full range 
of the frustration angle for different cluster sizes. Figure [T6l 
shows the maximum of the heat capacity as a function of ^. 
The bottom part of the figure shows the frustration dependence 
of the temperature Tcy at which the maximum is reached. 

Two overall effects are clearly visible: (1) Apart from the 
regime with strong frustration, the maximum increases with in- 
creasing cluster size. (2) The maximum temperature decreases 
with increasing cluster size. Taken together, this indicates that 
entropy is shifted to lower temperatures, a sign of the missing 
long-range coiTelations not included in the partition function 
for the small clusters. 

Our results are in qualitative agreement with those in fl, 
[Toll . They represent a quantitative improvement over the es- 
timates of f5]. Direct comparison with |19] is made difficult 
by the ambiguities associated with analytic continuation of a 
series using Pade approximants, and by the fact that the lim- 
ited number of cluster sizes we can use at present do not per- 
mit a finite size scaling analysis. In agreement with 01, 
drops sharply near the crossover between the spin liquid regime 
and the collinear phase around 72/7i w 0.6, corresponding to 
0/7r « 0. 17. Similar drops occur at the borders of the FM regime 
with the NAF and CAF phases, respectively. These drops are 
accompanied with a smaller Tcy in order to conserve the en- 
tropy of the system. 

In Figure the behaviour of the maximum of the mag- 
netic susceptibility J™" together with the temperature at which 
the maximum is reached is displayed. In contrast to the heat 
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Fig. 16. Maximum of the heat capacity Cv{T) and its position Tq, as 
functions of the frustration angle (j). The open (solid) circles denote 
the results for the 20-site (16-site) cluster, the dotted line denotes the 
eight-site cluster. 



Fig. 17. Maximum of the uniform magnetic susceptibility %{T) and 
its position as functions of the frustration angle <p. The open (solid) 
circles denote the results for the 20-site (16-site) cluster, the dotted 
line denotes the eight-site cluster. 
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capacity, X {T) does not display an anomaly upon crossing the Table 4. List of momentum vectors q for the 16-site square together 
spin-liquid regime. The maximum value diverges near the crossovejp'ith their multiplicity m. 
to the FM regime, while its position approaches T =0, which 
is the expected behaviour. Apart from that, the parameter de- 
pendence of the maximum position is qualitatively the same 
as for Tcy . 



Due to the sharp drop of the maximum of the heat capac- 
ity, the ratio of the two temperatures and Tcy shows a pro- 
nounced anomaly which has a strong dependence on cluster 
size in the spin liquid regime at ^/ttwO.IV {J2/J\ ~ 0.6). 
In Figure we have plotted Txl^Cy a function of ^. For 
0.3 < <^/k < 0.85, i.e., in the collinear phase, T^^^jTcy does 
not depend on 0. For negative values —1/2 < ^1% < —0.15 
(in the Neel phase), T^/Tcy depends roughly linearly on the 
frustration angle, providing an additional means to determine 
this angle uniquely from thermodynamic measurements alone. 
(Unfortunately, none of the three compounds discussed in this 
paper fall into this category.) 



3.3 Spin structure factor 

The finite temperature Lanczos approach also permits the direct 
evaluation of correlation functions. We consider here the static 
spin structure factor given by 



N 



,iq(R,-R,) 



(s,-s,). 



(37) 



We have calculated the temperature dependence of the spin- 
spin correlation functions (S,Sj) for the 16-site cluster and 
performed the necessary summations to determine ^(q, T). The 
irreducible triangle of the Brillouin zone of the 16-site cluster 
contains six points which are listed together with their multi- 
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Fig. 18. Ratio 7^/7cv of the maximum positions of the uniform mag- 
netic susceptibility x{T) and the heat capacity Cv{T) as a function 
of the frustration angle (j>. The solid (dashed) lines denote the results 
for the 20-site (16-site) cluster, the dotted line denotes the eight-site 
cluster. 

plicity in table |4] All points lie on the edges and corners of the 
triangle. 

In Figure^] S{q) is displayed as a function of q lying on 
the edge of the irreducible Brillouin zone triangle. We have 
chosen Ji/k^ = 1.25K, = 5.95K (collinear phase) and 

Ji/ks = 5.95 K, J2/kB 1.25K (Neel phase). These values for 
Ji and J2 correspond to those found for Li2VOSi04 in |19]; 
the former two coiTespond to (p^/n = 0.43, the latter have a 
frustration angle of (p+/n = 0.07. The dots represent the nu- 
merical results; the lines connecting them are just guides to the 
eye. We plot S{q,T) for ten different temperatures k^T/Jc = 
1 , 2, 3 ... 10 with an offset of one half between each two curves. 

For Pb2VO(P04)2, currently no suitable single crystals are 
available to be able to measure S{q, T) by diffuse neutron scat- 
tering. Therefore, we have calculated the angular average over 
the momentum transfer of S{q,T), which can be experimen- 
tally determined using powder or polycrystalline material. The 
results are displayed in Figure l20l as 5(|q|,r) versus the mod- 
ulus q = \q\ of the momentum transfer and the temperature T. 
We have chosen the two frustration angles 0_ = 0.64 where the 
system is in the collinear phase and 0+ = —0. 1 1 corresponding 
to the Neel phase. For the former, the maximum of 5'( |q| , T) is 
located at |q| = n, the latter reaches its maximum near the zone 
boundary where |q| = VlTt. 

To summarize, from Figures EH and m we can conclude 
that the structure factor provides a means to determine the frus- 
tration ratio (p and therefore the ordering wave vector unam- 
biguously: For S{q) is strongly peaked at q = (?r,0) or 
(0,7r), while for 0+ it is peaked at q = {n,n) as one would 
expect from the associated broken symmetries for in the 
thermodynamic limit. 

Figure [^displays the temperature dependence of S{q,T) 
for different values of q. In the collinear and the Neel phase, at 
low temperatures, S{q,T) develops a pronounced anomaly at 
the ordering vector Q characterizing the phase. This anomaly 
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Fig. 19. Static spin structure factor S{q,T) of the 16-site cluster for 
Jl/ke = 1.25 K, Ji/k^ ~ 5.95 K (collinear phase, top) and J\/k^ = 
5.95 K, Ji/k^ = 1.25K (Neel phase, bottom figure). The values cho- 
sen for Ji and Ji correspond to those determined for Li2VOSi04 
in 1 19]. The lines are guides to the eye, the dots denote the numeri- 
cal results. The individual curves in each figure, from bottom to top at 
q = F, correspond to fixed temperatures k^T jj^ = 1, 2, 4, and 8. 



is the precursor of a divergence of S{q, T ^Q) for the infinite 
system. The asymptotic value for k^T / ^ °° is S{q) ~ S{S + 
1) = 3/4 in each case, as it should be. 



In the collinear as well as in the Neel phase, the relation 
S{Q) > S{q), q 7^ Q holds for all temperatures, while for the 
spin liquid regime, this is qualitatively different: Here, the value 
of S{q,T) is approximately the same for q = {7t,n) and q = 
(tTjO) or (0,71:) at temperatures k^T > Jc and always larger 
than the value for q = 0, which is an additional indicator for 
the strong frustration in that phase. The approximate equaUty 
of S{q, T) values for different wave vectors also supports the 
picture that close to the SL phase domains of the CAF phase 
may easily form inside the NAF and vice versa as illustrated in 
Figure|3] 



To illustrate our numerical findings further, we have taken 
the Fourier transform of the high-temperature series expansion 
for the structure factor up to first order using the result in 
page 709ff] 
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Fig. 20. Static spin structure factor 5(|q|,r) of the 16-site cluster for 
-/l Ab = — 6K, ^2/^6 = lOK (collinear phase, top) and J\/k^ = lOK, 
Jl/ks ~ — 6K (Neel phase, bottom figure). The values chosen for Ji 
and J2 correspond to those given for Pb2VO(P04)2 in 0]. 
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(39) 



where z = 4 is the coordination number of the lattice. In Fig- 
ure 12] the dash-dotted lines in the three panels for the dif- 
ferent regimes in the (7i ,72) phase diagram represent the first 
two terms of equation (I38> for the corresponding characteris- 
tic wave vectors q = (7r,0) or (0,;r) in the collinear phase, 
q = {7i,n) in the Neel phase, and both of these values in the 
spin liquid regime. 

In the two phases where the system has a magnetically 
ordered ground state, the high-temperature approximation for 
5(q) at the respective ordering vector akeady underestimates 
the exact-diagonalization results for temperatures k^T/Jc w 5 



Fig. 21. Static spin structure factor 5(q, T) as a function of tempera- 
ture for q = (0,0) (solid line), q = (7r,0) or (0, k) (dotted line), and 
q = {n, n) (dashed line) for three different frustration angles (p = 7t/2 
(collinear phase, top), <j) = tan^'(— 1/2) (Neel phase, middle), and 
<j> — tan^'(l/2) (spin liquid, SL, phase, bottom figure). The dash- 
dotted lines in each of the three panels denote the high-temperature 
expansions (HTSE) of the structure factor at the respective value for 
q where S{({.T 0) diverges, i.e., from top to bottom q = (7r,0) 
for the collinear antiferromagnet, q = (n, n) for the Neel phase, and 
q = {7t,0) or {n, k) in the spin liquid regime. 



and below. In contrast, for the spin liquid regime, 5(q) is well 
reproduced for temperatures k^T well below J^, demonstrating 
that long-range coiTelation effects are suppressed due to the 
presence of the strong frustration. From equation ( I38> we can 
also conclude that at high temperatures kgT Jq the relation 
S{q = (;r,0)) = S{q = (;r,;/r)) holds exactly for J2/J1 = 1/2. 
We have chosen this particular frustration ratio for the lower 
panel in Figure |2] 
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Fig. 22. Angular average S{q, (j)) of the first-order term of the high- 
temperature series expansion of the static spin structure factor 5(q, T) 
as a function of the modulus q of the momentum transfer and the 
frustration angle (j) = —7t/2...K. (j) = —k/2 corresponds to Ji = 0, 
Jl/Jc = — 1> which is the border between the ferromagnetic and the 
Neel phase in FigureQ <j) = k corresponds to Ji //c = — 1 and J2 = 0. 
Going from (j) = —7t/2to (j) = 7t,we successively scan the Neel phase, 
the spin liquid regime, the collinear and finally the ferromagnetic 
phase (see text). 
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Fig. 23. Same as Figure l22l here for the two fixed values tj) = tj)± for 
Pb2VO(P04)2 (Table|3j. Solid line: 0/;/r = = 0.64 (collinear 

phase), dash-dotted line: (j)/n = (p+fn = —0. 1 1 (Neel phase). 

In Figures |22] and |23l we have plotted the angular average 
S{q,(p) of the first-order term of the high-temperature series 
expansion, equation i39\ . given by 

1 

S{q,^)^ — / da5i(q,0), q = ^(coso:,sina). (40) 
27t Jo 

Figure |22] demonstrates how the maximum of this coefficient 
of the structure factor evolves as a function of (j) . 

Figure l23lshows S{q,(p) as a function of the modulus of the 
momentum transfer q — |q| for the two possible values (j)± of 
the frustration angle for Pb2VO(P04)2. The maximum for 0+ 
is near the zone boundary, whereas the maximum for 0_ oc- 
curs at smaller values for q, confirming the conclusions of our 



exact-diagonalization results. Hence, for k^T ^ J^, it is pos- 
sible to determine experimentally the correct value of ^ from 
the quantity S{q,(p) ^ k^T /J^ {S{q, T) — S{S+ 1)), and there- 
fore together with the temperature dependences of xiT) and 
Cv{T) the exchange parameters Ji and J2- 

4 Summary and Conclusions 

Motivated by the discovery of Pb2VO(P04)2, a "J1-J2 com- 
pound" with at least one FM J, we have extended the semi- 
classical description of the J1-J2 model to the case with FM 
couplings. We discussed the possible nature of the phase tran- 
sitions between the three dominant phases, FM, NAF and CAP. 
On the basis of our results, the transition from CAF to FM 
seems to have much in common with the transition from CAF 
to NAF, where an intermediate spin liquid region is known to 
occur 

In addition to the phase diagram, the finite temperature prop- 
erties of the model were discussed in the light of numerical re- 
sults for 16 and 20 site clusters, using a newly implemented fi- 
nite temperature Lanczos algorithm. In particular, we used this 
to calculate spin-spin correlation functions at high temperature, 
which can be compared with diffuse neutron scattering experi- 
ments. Our numerical results should be of use in resolving the 
controversy which surrounds parametrising J1-J2 compounds. 

We also constructed a simple "tetragonal" mean field the- 
ory for the J\-J2 model, and performed an exact analytic di- 
agonalization of an 8-site cluster While these cannot be relied 
upon for quantitative comparison with experiment, they can be 
used to fit susceptibility data very easily, and seem to capture 
the essential physics of the model. Both are discussed in the 
Appendix. 

Much of the most interesting magnetic physics occurs in 
frustrated FM's - He III on graphoil and the CMR mangan- 
ites, to name but two - and it is our belief that the J1-J2 model 
deserves further study as such. 
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Fig. 24. The geometric content of the tetragonal mean-field theory for 
the 71 -J2 model. 
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Table 5. The energy spectrum of a single tetrahedron. 



H. Nishimori and Y. Ozeki, Journal of the Physical Society 
of Japan 58, 1027(1989). 



A Tetragonal Mean field theory 

We take as a starting point the Hamiltonian in equation (|8} for 
a single tetrahedron (cross-linked square), and treat this fully 
connected 4-site cluster as a building block for the lattice. From 
knowledge of the spectrum and degeneracies given in Table |5] 
we can can calculate the partition function of a single tetrahe- 
dron exactly, and from that, its magnetic susceptibility and heat 
capacity |33]. 

Starting from knowledge of the exact susceptibility of a 
tetrahedron we can construct a mean field theory for 

the lattice 



l+3(7i+72)Z™T(r) 



(41) 



The geometric content of the mean-field approximation is ex- 
plained in Figure |24] - one quarter of the tetrahedra making 
up lattice are treated exactly (double counting J\ bonds), and 
the remaining bonds are treated on a mean-field level in such a 
way that each spin sees three exact and nine mean field bonds. 
Because of the very high degree of frustration involved, this 
approach should provide a reasonable account of the high tem- 
perature susceptibility when the ground state of the tetrahedron 
is a singlet, i. e. in the bulk NAF and CAE phases. 

At high temperatures it is possible to expand the magnetic 
susceptibility found this way in \/T, and to make a direct com- 
parison with high temperature series expansions. Only the first 
two terms of equation (I32t are reproduced exactly, but the func- 
tional form of x{T) seems none the less to give a reasonably 
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Fig. 25. The enumeration of sites used for the 8-site cluster. 
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Table 6. The energy spectrum of the 8 site cluster 



good account of the experimental measured susceptibility. Since 
the partition function can be calculated for arbitrary magnetic 
field h, it is also easy to examine non-linear effects in x{T,h). 
The 4-site cluster is however, too small to give a convincing 
description of cy (T). 



B Exact solution of the 8 site cluster 





Fig. 26. The 9-site cluster with open boundary conditions used for the 
calculation of the lower bounds. 



where by ^2345678 denote the Heisenberg Hamiltonian on 
a complete graph spanned by the sites 1 to 8, i.e. 



^12345678 



E E S,S, = i i:S,- -^ES^ (43) 



= l7=i+l 



and similarly for the others. The nice feature of the Hamilto- 
nian (I42> is a particular hierarchy of the terms: starting from 
the basic building blocks - the states on sites (1,2), (3,4), 
(5,6), and (7,8), chosen either as singlet or triplet - we can 
construct successively the definite spin states on sites 1, 2, 3, 
and 4 (Suii) and 5, 6, 7, and 8 (55678), and finally on the whole 
cluster (5i2345678), and the energy for that state can readily be 
determined. So the problem is reduced to bookkeeping. The 
spectrum is shown in Table|6l 

^From the knowledge of the spectrum and degeneracies, 
we can calculate the exact partition function of the 8-site clus- 
ter, and from that its heat capacity and magnetic susceptibility, 
just as for the 4-site cluster. Expanding the magnetic suscepti- 
bility obtained in this way at high temperatures, we again find 
that we reproduce correctly only the first two terms of the high 
temperature series expansion. However empirically, the mag- 
netic susceptibility of the 8-site cluster provides an excellent fit 
to experiment. The estimate of heat capacity obtained in this 
way is also more reliable than that found from the 4-site cluster 
of SectionlAl 



In Section |X| we used the property that the energy spectrum 
of a complete graph depends only on the total spin to solve 
a fully connected 4-site cluster We can use a generalization 
of the same trick |34] to solve the 8-site cluster with periodic 
boundary conditions's shown in Figure|25l Written in terms of 
complete graphs, the Hamiltonian of this cluster is: 

-^sites = J 1 -^2345678 + (2-/2 - -^1 ) (-^1234 + -^678 ) 

-272 (^2° + ^34° + + ^78°) (42) 



C Lower-bound calculation 

A lower-bound estimate of the ground-state energy for spin sys- 
tems can be obtained by exactly diagonalizing the Hamiltonian 
on a small cluster with open boundary conditions. The small 
cluster is chosen such that it can be used to tile the complete 
lattice, in which case we can write the Hamiltonian of the sys- 
tem as 

M' = Y^Ml (44) 
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where is the Hamiltonian for an individual cluster. The 
simplest such cluster for the square lattice with nearest neigh- 
bour bonds is the cross-linked square (tetrahedron), discussed 
above. The ground state wave function of the complete lattice 

Jf|vAo)=£o|v^o) (45) 

can then be used as a variational wave function, providing an 
upper bound on the exact ground state energy iiciuster of a single 
cluster 

(VAo|=^|ro> >£cluster (46) 

It follows that ^cluster must in turn provide a lower bound on the 
ground state energy Eq per spin of the complete lattice. This ap- 
proach was introduced by Anderson |35], and later improved 
by Nishimori and Ozeki, who optimized the lower bound by 
treating the exchange on the bonds within the cluster as a varia- 
tional parameter 1 36]. We have refined the method further by al- 
lowing additional longer range bonds in the cluster which can- 
cel when the lattice is tiled by the sum over translated clusters. 
The 9-site cluster we used is shown in Fig.|26l with constraints 
AKi + 2K[ = 7i , 2K2 + 2K'2 = J2 and 2Ko, +K'-^=Q. 



